Re-integerization of fractional charges in the correlated quarter-filled band 



O 

o 
o 

(N 

> 

o 
O 

oo 



i 

CO 

-4-5 
c3 



I 

C 

O 

o 



> 

(N 
l> 
(N 
O 

o 
o 

c3 



c 

o 
o 



X 
S3 



R.T. Clay 1 ' 2 , S. Mazumdar 1 , and D.K. Campbell 3 
1 Department of Physics, University of Arizona, Tucson, AZ 85721 

2 Cooperative Excitation Project ERATO, Japan Science and Technology Corporation (JST), University of Arizona, Tucson, 

AZ 85721 

3 College of Engineering, Departments of Electrical and Computer Engineering and Physics, Boston University, Boston, MA 

02215 
(February 1, 2008) 

Previous work has demonstrated the existence of soliton defect states with charges ± e/2 in the 
limits of zero and of infinite on-site Coulomb interactions in the one- dimensional (ID) quarter-filled 
band. For large but finite on-site Coulomb interaction, the low temperature 2kF bond distortion 
that occurs within the 4k f bond-distorted phase is accompanied by charge-ordering on the sites. We 
show that a "re-integerization" of the defect charge occurs in this bond-charge density wave (BCDW) 
state due to a "binding" of the fractional charges. We indicate briefly possible implications of this 
result for mechanisms of organic superconductivity. 



PACS indices: 71.30.+h, 71.38.+i, 71.45.Lr, 74.70 
One of the most dramatic and exciting predictions 
of the initial model studies of quasi-lD organic 7r- 
conjugated polymers and charge transfer solids (CTS) 
was the existence of exotic excitations, including solitons 
with unusual spin-charge relations and fractional 

charge For trans-polyacetylene, with a 1/2-filled 

ID band, it was shown within the Su-Schricffcr-Heeger 
(SSH) one-electron model that the neutral solitons have 
spin 1/2, while charged solitons are spinless |l|,||. This 
result remains unaltered upon inclusion of the Coulomb 
electron-electron (e-e) interactions [9|JlO|]. In the case of 
the ID 1/3- and 1/4-filled bands, any non-zero electron- 
phonon (e-ph) coupling leads to 2kp Peierls transitions 
corresponding to trimerization and tetramerization of the 
respective lattices, and defect soliton charges are ± 2e/3 
HQ and ± e/2 [||, respectively. Effects of e-e interac- 
tions on the soliton charge in non-l/2-filled bands have 
been investigated only for the 1/4-filled band, primarily 
for the limiting case of large on-site Hubbard repulsion 
(U -> oo) [| §. For U -> oo, the ID 1/4-filled Hubbard 
band of electrons reduces to a 1/2-filled band of nonin- 
teracting spinless fermions . Hence in the presence of 
an e-ph coupling, a Peierls transition occurs at twice the 
fermi wave vector corresponding to the spinless fermions, 
which is four times the original kp (kp = 7r/4a, where 
a is the undistorted lattice constant). Depending on the 
nature of the phonons, this leads to bond or site dimer- 
ization of the lattice. A single doped hole/electron gener- 
ates two charged soliton defects in the dimerized lattice, 
which therefore have charge ± e/2 again. Since soliton 
charges are e/2 at both [7 = and U — > oo, a continuity 
between these two limits has been claimed ||. 

In the case where the Akp metal-insulator transition is 
due to bond distortion, a 2kp spin-Peierls (SP) dimeriza- 
tion of the dimerization (i.e., tetramerization) occurs at 
low temperatures for finite U [[llj . Reference || , which 
claimed the persistence of fractional charge in the 2kp 
phase, did not actually probe the nature of the 2kp. We 
have recently shown that the 2kp SP bond tetrameriza- 
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tion is also accompanied by an on-site 2kp CDW |l2"| , 
and the correct description of this state is a bond-charge- 
density wave (BCDW). Clearly an important theoretical 
question is whether fractionally charged solitons continue 
to exist in the BCDW ground state. 

This issue is of considerably more than mere theoret- 
ical interest in the case of the organic CTS. The origi- 
nal applications of the "large-C/" 1/4-filled band theory 
were to 1:2 anionic TCNQ solids fpL3f| , where, however, 
no direct evidence for fractional charge was found, in 
spite of experimental attempts. More recently, interest 
has shifted to the nominally 1/4-filled 2:1 cationic CTS, 
many of which exhibit superconductivity and other novel 
phases at low temperatures [Q. The most strongly ID 
systems here, (TMTTF) 2 X (X = PF 6 , AsF 6 etc.) ex- 
hibit a high temperature metal-insulator (presumably 
4kp) transition, and a SP transition at T$p < 20 K. 
Again, no evidence for fractionally charged defect states 
has emerged. The question that arises then is whether 
the half-integer charges are unobservable, or the right 
experiments have not been done yet. This question ac- 
quires added importance when the entire family of 2:1 
cationic CTS is considered. In our recent theoretical 
studies, for example, we have shown that the BCDW 
persists in the interacting 1/4-filled band for all t±_ < t, 
where t and t± are the intra- and interchain one-electron 
hopping integrals, respectively, and furthermore, in the 
small t± region the BCDW also coexists with a spin- 
density wave (SDW), giving a BCSDW (l|Jl|. Ex- 
perimentally, charge-ordering [l^Jl^| has been found in 
the ID (TMTTF) 2 X which are SP systems, while the 
BCSDW has been seen in (TMTSF) 2 X f^fo) and a- 
(BEDT-TTF) 2 MHg(SCN) 4 @. Evidence for coexisting 
SP-like states and charge-ordering is also beginning to 
emerge in other BEDT-TTF systems @^|. Based on 
the persistent BCDW that is obtained theoretically and 
also observed experimentally, we have suggested that su- 
perconductivity in these strongly correlated systems is 
a consequence of pairing of commensurability defects of 
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the BCDW phase |T(| . If these commensurability defects 
are fractionally charged then, an obvious question that 
emerges is whether the the superconducting "pairs" are 
actually "quadruplets" ! 

In the present Letter, we show that charge- 
fractionalization does not occur in the low temperature 
2kp SP phase of the interacting ID 1/4-filled band: the 
charge-ordering associated with the SP tetramerization 
causes a "re-integerization" of the fractionally charged 
defects that occur within the 4kp dimerized phase. This 
re-integerization provides support to the idea that super- 
conductivity in the quasi-2D TMTSF and BEDT-TTF 
may be due to pairing of integer charged commensura- 
bility defects in a background BCDW, although it does 
not necessarily exclude the "melting" of the BCDW at 
the superconducting transition. 

To establish this result, we will consider the ID Peierls- 
extended Hubbard Hamiltonian for the case of bandfilling 
slightly away from 1/4, and will present both many-body 
numerical results as well as heuristic physical arguments. 
No SP transition can occur within the 4kp CDW [fL6| , 
and therefore the the experimentally observed SP transi- 
tion in (TMTTF) 2 X proves that the 4kp metal- insulator 
transition is to a bond-order wave (BOW) that occurs 
for intersite e-e interaction V < V c ~ 2\t\ p6[ . Upon 
the exclusion of the 4kp charge-ordering, there are four 
possible equivalent charge-orderings that can occur for 
nonzero U and V < V c , viz., ...1100..., ...0110..., ...0011..., 
and ...1001... (here and below a '1' implies charge 0.5 + 
e, and a '0' implies 0.5 - e). The 4k^ bond dimerization 
breaks a twofold symmetry and the ground state is a su- 
perposition of the type ...1100... + ...0011..., with weak 
(W) '11' and '00' bonds, and strong (S) TO' and '01' 
bonds. The 2kp SP transition breaks yet another twofold 
symmetry, and the ground state is now either ...1100... 
or ...0011..., with the TP (W) and '00' (W) bonds now 
inequivalent, and the bond distortion pattern W'SWS 
|P^| , [l5|Jl6| . Our goal is to calculate defect charges within 
this W'SWS state. 

Since our many-body numerical results are for finite 
chains, it is useful to establish our methodology by first 
solving the known case of U — > oo, V = 0. The spinlcss 
fermion Hamiltonian here is written as, 

H SF = -t^2(a]a j+1 + h.c) (1) 
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For simplicity, we have not included the SSH e-ph in- 
teraction that gives the bond dimerization in the above 
HI). As we show below, calculations of the bond orders, 
defined as (Bjj + i) = (ajaj+i + h.c), and of site charge 
densities, of an open long finite chain with uniform hop- 
ping t show spontaneous modulations of the bond orders 
and charge densities, with the same modulation patterns 
as expected within a fully self-consistent calculation of a 
periodic electron-phonon coupled ring. The actual lat- 



tice dimerization occurs for + electron-phonon coupling 
in the infinite chain. 

Figure [t](a) shows the bond orders of a 64-site open 
chain of noninteracting spinless fermions with uniform 
hopping integrals; the spontaneous period-two dimeriza- 
tion of the chain is evident in spite of end effects. In 
Figure |l|(b), we show the result of adding two holes. The 
two holes produce four defects (regions where consecutive 
bond orders are equal), each of which then have charge 
+ e/2. Charge density calculations (not shown) give uni- 
form (rij) = 0.5 for the undoped chains, and four defects 
in the doped chain, with maximum deviations from 0.5 
ocurring exactly in the regions where consecutive bond 
orders are equal. 
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FIG. 1. The bond-orders of a 64-site uniform open chain of 
noninteracting spinless fermions with (a) 32 fermions, show- 
ing the spontaneous dimerization that occurs for the effective 
half-filled band; (b) two doped holes (30 fermions), showing 
the resulting four defects. 

Although the above figures reproduce known results, 
it is useful to give a configuration space interpretation 
of these results. In Figure |^(a) we show the sponta- 
neously dimerized ground state, with equal (rij) = 0.5 
on each site. Here the bonds between the sites within a 
box are stronger than the bonds between sites that oc- 
cur in neighboring boxes. Removal of a single electron 
initially removes two charges 0.5 from each of two adja- 
cent sites connected by a weak bond (see first diagram, 
Figure || (b)). The configuration with adjacent defect 
sites has higher energy than the one in which they are 
separated [jl],|j , which leads to a spontaneous separation 
of the defects (see second diagram, Figure | (b)), each 
carrying charge e/2. The defects are shown localized to 
a single site in Fig. |](b) for simplicity; as seen in the 
numerical data of Fig. |](b), each defect actually spreads 
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out over several sites. 




FIG. 2. (a) The undoped 4k F dimerized 1/4-filled band. 
Boxes represent the "dimers", i.e, sites coupled by strong 
bonds, (b) creation of two e/2 solitons upon the removal 
of one electron from the configuration of (a), (c) The BCDW 
state that describes the interacting 1/4-filled band, below the 
2k_p SP transition (see text), (d) A doped hole removes a 
single '1', and soliton charge is now e. (e) Configuration with 
two doped holes in the 2ki? BCDW. 

To explore the defect charges for finite U below the 
2kp transition, we study the Hamiltonian Pq,E3[: 



if = -*1 2J B 2 j-l,2j,a - X^[ <2 ~~ a/ ^j} B 2j,2j+l,cr (2) 
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In the above B%^,a = c\ a Cj. a + h.c. is the standard hop- 
ping operator for an electron with spin a, nj >(7 — cj a c^ a 
and rij — ^ CT the standard electron-number oper- 
ators, a and K the electron-phonon coupling and the 
spring constants, respectively, and Aj is the displace- 
ment of each dimer from equilibrium. The intra-dimer 
hopping t\ is larger than the inter-dimer hopping t 2 be- 
cause of the 4ki? bond dimerization which is assumed to 
have occurred already. 

Before showing our numerical data, we present physi- 
cal arguments based on a configuration space picture that 
already suggest re-integerization of the fractional charge 
within the above Hamiltonian. As shown in our previous 
work (see also below) , the ground state of the commensu- 
rate interacting 1/4-filled band is the BCDW state shown 
in Figure |](c). Here each TO' and '01' within a box are 
connected by the strong (S) bond, while the "double" 
and "single" bonds correspond to the W and W bonds, 
respectively. The ...1100... charge-ordering dictates that 
"doping" now will involve removal of a whole electron 
from a single site, as shown in the first diagram of Fig- 
ure |^(d). Simply removing the T' creates an energetically 
unfavorable asymmetric situation with three '0's in a row, 
and energy is gained by the reorganization shown in the 
second step of Figure |](d), which reduces the number of 
weak '0-0' bonds. The difference from the U — ► 00 case 
is then as follows. The defect induced charge deficiency 
now occupies a single box and no further separation be- 
tween the charges can occur. Clearly, the soliton charge 
is then +e and not +e/2. Note that the charge reorgani- 
zation creates phases of inter-dimer bond dimerizations 



that are different on the left and the right of the soliton 
defect (this may be seen most easily by numbering the 
boxes sequentially). Soliton defects with integer charge 
therefore occur in pairs, where a second hole "heals" the 
"phase-problem" induced by the first hole, as shown in 
Figure |](e). Note that the above physical picture, though 
qualitative, predicts two important numerically testable 
details, viz., (i) the intra-dimer bonds become weak near 
the defect center, with a minimum in the intra-dimer 
bond order occurring over a single dimer, (ii) the inter- 
dimer bonds around a defect increase in strength, and 
reach their maximum values on both sides of the weakest 
intra-dimer bond. 

In our numerical investigation we study the bond or- 
ders J2a(Bj,j+i,<r) an d charge density X) CT ( n J>)' a * an< ^ 
near 1/4-fhling. We use the Stochastic Series Expan- 
sion (SSE) quantum Monte Carlo (QMC) method [§4|, 
which is exact to within statistical error for this ID sys- 
tem: there is no approximation in the discretization of 
imaginary time, as in many other QMC methods. As 
SSE is a finite-temperature QMC method, a sufficiently 
low temperature must be used to obtain ground-state re- 
sults. In the results shown below, inverse temperatures 
(3 = 4L (L=chain length) were used, which showed es- 
sentially no difference from (3 = 2L computations. Since 
bond-distortions in open chains decrease with length, un- 
like in periodic rings, where they increase (convergence to 
the same finite amplitude for the BOW occurs in the long 
chain limit) our initial calculations arc for relatively short 
chains with 32 sites. Figure ||(a) shows the bond orders 
for a 32-site open chain with 16 electrons and U = 6, 
V = 1, and hopping integrals t\ — 1.4 and t 2 = 0.6 
(in units of t — (ti + t 2 )/2). It is seen that the strong 
intra-dimer bonds have nearly uniform bond-orders, but 
there occurs a spontaneous alternation of the inter-dimer 
bond orders, giving an overall bond distortion pattern 
W'SWS. Figure ||(b) shows the charge densities in this 
undoped 1 /4- filled system, and there is clear spontaneous 
site-CDW with periodicity ...1100... in the center of the 
chain, with relatively weak chain end effects. These re- 
sults show that once the 4k^ bond dimerization (t± ^ t 2 ) 
is explicitly put in "by hand" , further dimerization of the 
dimerization is unconditional. 

In Figure ||(c) we show the bond orders for the dimer 
chain "doped" with two holes. There occur now only 
two defects, centered on the dimer bonds between sites 9 
and 10, and sites 23 and 24, respectively. These defects 
therefore have integer charges. Note that the intradimer 
bond orders are smallest for the 9-10 and 23-24 bonds, 
while the interdimer "weak" bonds on both sides of the 
defect centers reach their strongest values, in complete 
agreement with Figures ^(d) and (e). 

The very large difference between t\ and t 2 in Fig- 
ures |3|(a) - (c) was chosen merely to make the results 
visually most accessible, and is not a required condition 
for obtaining defects with integer charge. To demonstrate 
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FIG. 3. SSE bond orders of 32 site open chains with {7 = 6, 
V = 1. The lines are guides to the eye. (a) Undoped dimer- 
ized 1/4-filled system, ii = 1.4 and t 2 — 0.6. Spontaneous 
dimerization of the dimer lattice results in the alternation 
of weaker bond orders, (b) Charge densities of the undoped 
chain, showing ...1100... ordering, (c) Bond orders with two 
added holes, with the same parameters as in (a), (d) Same 
as (c), but with weaker intrinsic dimerization, ti = 1.2 and 
t 2 = 0.8. 



this, we have performed calculations with smaller ti—t%. 
In Figure ||(d) we show our results for the same 32-site 
system with two doped holes, but now for weaker 4k^ 
dimerization t\ — 1.2 and t 2 — 0.8. No change in the 
defect structure is seen with the weaker dimerization. 

Finally, it is important to demonstrate that the occur- 
rence of only two defects in Figures ||(c) and (d) is not 
a consequence of short chain lengths and overlapping de- 
fects. We have therefore performed our calculations also 
for 64 sites, a chain length where C/ — > oo shows clear ev- 
idence for four defects (see Figure ||) . Based on the simi- 
larity between large and moderate t\ — 1 2 in Figures. ||(c) 
and (d) , we performed our calculations here for only large 
ti — *2j viz-, ti = 1-4, t 2 — 0.6. Our SSE results for a 
64-site open chain with 30 spinful electrons are shown in 
Figure ^, where once again only two defects (centered on 
the dimer bonds between sites 17 and 18, and sites 45 
and 46, respectively |2^]) are seen. An important feature 
that is clearer in this longer chain is that the weaker bond 
orders correctly recover their strong- weak character away 
from the defect centers. 

To conclude, the apparent continuity of fractional 
charges between the U = and U — > oo limits in the 1/4- 
filled band p|] arises only when the existence of the 2ki? 
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FIG. 4. SSE bond orders of a 64 site chain with two added 
holes, for U = 6, V = 1, ti = 1.4, t 2 = 0.6. 



BCDW ground state at finite U is ignored. The charge 
e/2 solitons that occur at U = are domain walls be- 
tween four distinct density wave phases, while at U — > oo 
they are domain walls between only two such phases. 
The e/2 charge at U — > oo is thus a consequence of site- 
occupancies of 0.5 electrons alone, and re-integerization 
of the fractional charge occurs in the 2kj? BCDW ground 
state. The integer charge solitons do not bind in ID, 
but it is conceivable that integer charge commensurabil- 
ity defects in the 2D BCDW bind into large bipolarons, 
thereby giving the superconductivity observed in the 2D 
CTS. Recent demonstrations of large bipolarons [^6| with 
strongly interacting 2D electrons, and the finite mobility 
of such intersite bipolarons 1 26 27j are highly promising 
results in these context. 

RTC acknowledges helpful discussions with AW. 
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